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Abstract 

Dynamic phase transition properties of ferromagnetic thin film system under the 
influence both bias and time dependent magnetic fields have been elucidated by 
means of kinetic Monte Carlo simulation with local spin update Metropolis 
algorithm. The obtained results after a detailed analysis suggest that bias field 
is the conjugate field to dynamic order parameter, and it also appears to define 
a phase line between two antiparallel dynamic ordered states depending on 
the considered system parameters. Moreover, the data presented in this study 
well qualitatively reproduce the recently published experimental findings where 
time dependent magnetic behavior of a uniaxial cobalt films is studied in the 
neighborhood of dynamic phase transition point. 

Keywords: Ferromagnetic thin film, dynamic phase transitions, Monte Carlo 
simulation. 


1. Introduction 

When a magnetically interacting spin system with ferromagnetic cou¬ 
pling J is exposed to a magnetic field oscillating in time, the system may not 
respond to the external perturbation simultaneously, and hereby, two impor¬ 
tant striking phenomena occur: Non-equilibrium phase transitions and dynamic 
hysteresis behavior. A typical ferromagnetic system exists in dynamically disor- 
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dered phase, in which the time dependent magnetization of the system oscillates 
around zero value in the range of high temperature and high amplitude of field. 
In contrary to the aforementioned treatment, for the low temperature and small 
amplitude of field regimes, the time dependent magnetization exhibits an oscil¬ 
lating behavior in a restricted range around one of its two non-zero values. For 
the first time, this type of investigation regarding the dynamic phase transi¬ 
tion properties of Ising model under the presence of a sinusoidally oscillating 
magnetic field have been performed by Tome and Oliveira by benefiting from 


Molecular Field Theory (MFT) [l|. Since then, a great deal works concerning 
the non-equilibrium phase transitions as well as hysteresis behaviors of different 
types of magnetic sjgtems have been investigated by a variety of techniques such 

23 . For ex- 


as MFT 


flIlflQ, 


and Monte Carlo simulations (MC) 


Effective Field Theory (EFT) 
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ample, thermal and magnetic phase transition properties of the kinetic Ising 


model have been analyzed within the framework of MFT in [6| , and it is found 
that frequency dispersion of the hysteresis loop area, the remanence and the 
coercivity have been categorized into three distinct types for varying system 
parameters. Furthermore, Park and Pleimling have considered kinetic Ising 
models with surfaces subjected to a periodic oscillating magnetic field to probe 
the role of surfaces at dynamic phase transitions. They have reported that the 
non-equilibrium surface universality class differs from that of the equilibrium 
system, although the same universality class prevails for the corresponding bulk 


systems 


fl¬ 


it is obvious from the above picture that only influences of a time depen¬ 
dent magnetic field on the dynamic phase transitions properties of the kinetic 
Ising model and its derivations have been realized, and possible situations have 
been addressed by means of theory and modeling. Besides, as far as we know, 
there are a rather few experimental and theoretical studies associated with the 
discussion of possible role of an additionally time-independent magnetic field, 
namely bias field hb, on the dynamic phase transitions in the magnetic systems 
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27|, and these studies suggest that hb appears to be conjugate 
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field of the dynamic order parameter. Very recently, this fact has been verified 
by experimentally and theoretically in Ref. 28| where time dependent mag¬ 


netic behavior of a uniaxial cobalt films under the presence of both bias and 
time dependent magnetic fields has been studied, and after a detailed analysis, 
it is observed that the bias field is conjugate field of the dynamic order param¬ 


eter. Berger and co-workers in Ref. 


28( 1 have used the MFT as a theoretical 


tool. It is a fact that spin fluctuations are ignored and the obtained results 
do not have any microscopic information details of system in MFT. Keeping 
in this mind, in order to qualitatively reproduce the experimental observations, 
we have implemented a series of MC simulations of a ferromagnetic thin film 
system described by a simple Ising Hamiltonian including both bias and time 
dependent oscillating magnetic fields. Some outstanding results are given in this 
letter, and it can be said that our findings qualitatively support and confirm 
the experimental results [^ . 

The paper organized as follows: In section [51 we briefly introduce our model 
and MC simulation procedure. The results and discussion are presented in 
section m and finally section [4] includes our conclusions. 


2. Formulation 

We consider a ferromagnetic thin film with thickness Lz under the exis¬ 
tence of both bias and time dependent magnetic fields. The Hamiltonian of the 
considered system can be written in the following form: 


H = -Y,,h,S,Sj-h{t)Y,S,, ( 1 ) 

where Si is conventional Ising spin variable which can take values of Si = ±1, 
and Jij is the nearest-neighbor spin-spin interaction term, and it is kept fixed as 
Jij = J(> 0) in throughout the simulation. The first summation in Eq. (|T]) is 
over the nearest-neighbor site pairs while the second one is over all lattice sites 
in the thin film system. h(t) denotes the time dependent oscillating magnetic 
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field term, and it composes of both bias and time dependent oscillating magnetic 
fields which has the following form: 

h{t) = hb + ho sm{ujt), (2) 


here hf, is the bias field, t is time, ho and w are amplitude and angular frequency 
of the driving magnetic field, respectively. The period of the oscillating magnetic 
field is given by t = 27r/a;. 

We simulate the system specified by the Hamiltonian in Eq. o on a 
Lx X Ly X Lz simple cubic lattice, and we apply free boundary condition in 
the z— direction of the thin film system, whereas in the directions perpendicu¬ 
lar to the z— direction we use periodic boundary conditions. We have studied 
the thin film system with thickness Lz = 20 with Lx = Ly = 100. The simu¬ 
lation procedure we follow in this work is as following: The simulation begins 
at a high temperature using a random initial condition, and then the system is 
slowly cooled down with a reduced temperature step fc^AT/J = 0.03, where 
ks and T are the Boltzmann constant and temperature, respectively. The con¬ 
figurations were generated by selecting the sites sequentially through the lattice 
and making single-spin-flip atternpts, which were accepted or rejected accord¬ 
ing to the Metropolis algorithm [29|, l30| . The numerical data were generated 
over 50 independent sample realizations by running the simulations for 20000 
MC steps per site after discarding the first 10000 MC steps. This amount of 
transient steps is found to be sufficient for thermalization for the whole range 
of the parameter sets. 

Our program calculates the instantaneous value of the total magnetization 
M{t) at time t as following: 


= NsMs{t)+NBMB{t) 

Ns+ Nb ’ ^ ^ 

where Ms{t) = and MB^t) = while Ns{= LxXLyX 2) 

and Nb{= Lx x Ly x {Lz — 2)) terms correspond to the total number of the 
spins located on the surface and bulk of the thin film system, respectively. In 
this study, we are interested only in thermal variation of total dynamic order 
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parameter as a function of system parameters, hence, from the instantaneous 
magnetization, we can obtain the dynamic order parameter as follows [l|: 



( 4 ) 


In order to specify the dynamic phase transition point at which dynamically 
ferromagnetic and paramagnetic phases separate from each other, we use and 
check the thermal variation of dynamic heat capacity which is defined according 
to the following equation: 

where E is the energy per spin over a full cycle of the external applied magnetic 
field which has the following form: 



( 6 ) 


3. Results and Discussion 

In this section, first of all, in order to understand the dynamic evolution 
of the magnetic system in detail, we will focus our attention on non-equilibrium 
phase transition properties of the ferromagnetic thin film system under the ex¬ 
istence of only a time dependent oscillating magnetic field. We will argue and 
discuss how the amplitude and period of the field affect the dynamic critical 
nature of the system. Next, for some selected combinations of Hamiltonian 
parameters, we will give and examine the bias field effects on the thermal and 
magnetic properties of the thin film system. Finally, we will discuss the com¬ 
peting mechanism between bias and time dependent magnetic fields. 

Dynamic phase diagram separating dynamically ordered and disordered phases 
of the thin him system in a {ho/ J — ksTc/J) plane in the absence of bias held 
is shown in Fig. [I] at various r values. At hrst sight, it can be easily seen 
from the hgure that, for a hxed value of r, the system shows ordered phase at 
the relatively low temperature and applied held amplitude regions. Here, the 
time dependent magnetization of the system oscillate around a non-zero value, 
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Figure 1: Dynamic phase diagram in a (/iq/J — kgTc/J) plane obtained from the peaks of 
the dynamic heat capacity for selected values of applied field period such as r = 50,100 and 
200 . 

namely, the system may not respond to the external magnetic field instanta¬ 
neously. With increasing strength of the amplitude of the external field, the 
magnetic phase of the system tends to shift dynamically paramagnetic phase. 
In this region, the time dependent magnetization of the thin film system oscil¬ 
late around a zero value, and it can follow the external field with a relatively 
small phase lag depending on the studied system parameters. One of the ob¬ 
tained results of our MC simulation study is that the thin film system shows 
sensitivity to varying applied field period in accordance with the expectations. 
For a fixed value of Hq/J, it is obvious that as r is increased, the phase tran¬ 
sition point is lowered because of the fact that decreasing field frequency leads 
to a decreasing phase delay between the magnetization and magnetic field and 
this makes the occurrence of the dynamic phase transition easy. The physical 
discussions mentioned above can be easily visualized by checking any two values 
of oscillation periods such as r = 50 and 200 for fixed value of ho/J = 1.5. 



6 





1.0 

B 

I 0.8 

ro 

L_ 

cc 

‘ 2 : 0.6 
0 
■Q 

o 

o 0.4 

E 

CO 

c 

Q 0.2 


0.0 

0 12 3 

ksT/J 


Figure 2: Effects of the external applied field period on the temperature versus dynamic order 
parameter curves corresponding the dynamic phase diagram illustrated in Fig. ^ 


In Fig. m we present the thermal variations of dynamic order parameter of 
thin film system corresponding to phase diagrams plotted in Fig. [T]for various r 
values with a selected applied field amplitude /iq/J = 1.5. It is evident from our 
simulation that when the temperature increases starting from zero, the dynamic 
order parameter starts to decrease from its saturation value, and it undergoes 
a second order dynamic phase transition between dynamically ordered and dis¬ 
ordered phases. As we discussed above, dynamic phase transition point moves 
downward in the temperature space with increasing applied field period. It is 
possible to see more clearly these facts in this figure. 

Before we start to discuss the influences of bias field on the thermal and 
magnetic behavior of the thin film system, it is beneficial to underline that there 
exists a competition between bias and time dependent magnetic fields. Namely, 
the bias field tries to keep the system in ordered phase {Q ^ 0) whereas the time 
dependent magnetic field tries to drive the system into dynamically disordered 
phase {Q = 0) depending on the other system parameters. From this point 
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Figure 3: For considered values of h^jJ = l.^^ksTjJ = 1.35 and r = 200, time series of 
magnetization of ferromagnetic thin film system, (a) and (b) are given for = 0.0 and 

0.15, respectively. 


of view, in Fig. [3] (a) we give time series of the magnetization without a bias 
field for selected values of h^/J = = 1.35 and r = 200. One can 

deduce from the figure that the unbiased magnetization curve displays a nearly 
rectangular shaped time sequence around zero value such that the thin film 
system exists in dynamically paramagnetic phase. As shown in Fig. Eljb), by 
keeping the system parameters fixed, if one applies the bias field to the thin 
film system, for example hf,/J = 0.15, the biased magnetization curve shows 
a small variation around its saturation value with increasing time. We should 
note that application of the bias field to the system is sufficient to induce a large 
Q value indicating the existence of an ordered phase, and it also allows us to 
analogy with static paramagnetic saturation at large applied field values. If one 
compares the our MC simulation findings with the experimental data in Fig. 2 


of Ref. 


■ 3 : 


a well agreement is seen between modeling and experiment. 


In the following analysis, we give the bias field dependencies of the dynamic 































Figure 4: Bias field dependencies of the dynamic order parameter obtained for values of 
k^T/J = 1.35 and h^jJ = 1.5. The dynamic hysteresis curves are calculated for different 
values of applied field periods: (a) r = 50, (b) r = 100 and (c) r = 200. The decreasing and 
increasing branches of the curves are denoted by downward and upward arrows, respectively. 


order parameter for a fixed temperature ksT/J = 1.35 and applied field am¬ 
plitude /iq/J = 1.5 for varying applied field periods in Fig. |4l[a)-(c). The mea¬ 
surement procedure to obtain the decreasing branches of the hysteresis curves 
is as follows: The bias field starts at -1-1.0/J, and then it is slowly decreased 
with a step Ah^/J = —0.02 until it reaches to —1.0/J. We obtain the increas¬ 
ing branches of the curves using a similar way. The decreasing and increasing 
branches of the hysteresis curves are shown by downward and upward arrows 
in figures. One of the outstanding results is that dynamically ordered state is 
not simply inverted at hj,/ J = 0.0, and therefore, Q — hb/J hysteresis curve ap¬ 
pears. When one compares the obtained curve with M — H hysteresis curve for 
ferromagnetic system at thermal equilibrium, it may possible to say that hb/J 
bias field appears to be the conjugate field of the dynamic order parameter Q. 
Furthermore, the aforementioned behavior strongly depends on the applied field 
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period, and it disappears with increasing value of the applied field period. 




Figure 5: For a fixed temperature value kgT/J = 1.35, (a) displays the stable dynamic 
order parameter for h^/J = 1.5 while (b-d) represent the variation of AQ in 2D contour 
representation in a (r—J) plane for varying applied field amplitudes such as (b) /iq/J = 1.5, 
(c) ho/J = 1.7 and (d) Kq/J = 1.9. 


As a final investigation, we touch upon the period dependencies oi ht/ J —Q 
dynamic hysteresis curve for a considered value of temperature fc^T/J = 1.35 
with different values of the applied field amplitudes in Fig. (I5|). In order to 
generate the Fig. [Sl[a) constructed as a color-coded map displaying the stable 
dynamic order parameter as a function of the applied field period and bias field. 


we follow a similar way used in Ref. 


It is found that a sharp transition 


takes place at = 0.0 between positive and negative Q values up to a 

particular value of r. With increasing value of r starting from a relatively 
low value, Q vanishes at h^/ J = 0.0. As a result of this, the sharp transition 
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behavior tends to disappear, and there is no phase transition line as t reaches 
at a critical value. The curves in Figs. [Sl[b)-(d) are plotted for three values of 
applied field amplitude: (b) ho/J = 1.5, (c) /iq/J = 1.7 and (d) h^/J = 1.9. 
The figures represent a map of AQ{hb/J,T) = Q‘^{hb/J,T) — Q^{hb/J,T). In 
other words, this formula determines the differences of the Q values obtained 
for the decreasing branch Q‘^{hb/J,T) and the increasing branch Q'‘{hb/J, t). 
It is clear from the figures that the obtained AQ{hb/J, t) curves appear as 
the roughly triangular structure extending up to a critical period value which 
sensitively depends on the applied field amplitude and other system parameters. 
When the energy coming from the time dependent applied field is increased, the 
boundaries of the triangular shapes tend to show a decreasing trend, and it is 
possible to mention that after a special value of the amplitude of field, the 
triangular shape completely disappears. 

4. Conclusion 

In summary, we have investigated the thermal and magnetic phase transi¬ 
tion features of the ferromagnetic thin film system under the influence both bias 
field and time-dependent magnetic fields. For this purpose, we have used Monte 
Carlo simulation technique with single-spin flip Metropolis algorithm. First, to 
determine the dynamic phase transition point, we have treated the thermal vari¬ 
ation of the heat capacity for considered values of the system parameters, and 
by benefiting from its peak, we have constructed the dynamic phase diagram 
in reduced magnetic field versus temperature plane. It is reported that for a 
fixed value of ho/J, as r is increased, the phase transition point is lowered since 
decreasing field frequency gives rise to a decreasing phase delay between the 
magnetization and magnetic field and this makes the occurrence of the dynamic 
phase transition easy. Next, we have focused on effects of the varying bias field 
on the ferromagnetic thin film system, and we have studied the dynamic hys¬ 
teresis behaviors as a functions of system parameters in detail. It is emphasized 
that hb is the conjugate field to the dynamic order parameter, and also hb ap- 


II 


pears to define a phase line between two antiparallel dynamic ordered states. As 
a final note, it should be underlined that our MC simulation results obtained in 
the present work corroborate the experimental observations found in Ref. H- 
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